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Abstract 



A statistical model is advanced for describing quantum turbulence in a superfluid 
C^ ■ system with Bose-Einstein condensate. Such a turbulent superfluid can be realized for 

trapped Bose atoms subject to either an alternating trapping potential or to an alter- 
nating magnetic field modulating the atomic scattering length by means of Feshbach 
^ ■ resonance. The turbulent system is represented as a continuous mixture of states each 

O . of which is characterized by its own vorticity corresponding to a particular vortex. 
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1 Quantum turbulence 

Turbulence in classical fluids has been intensively studied for many years [1-5]. It can also 
occur in superfiuids, such as ^He, ^He, and trapped atomic gases. A superfluid differs from 
a classical fluid in three principal respects: it exhibits two-fluid behavior; the superfluid 
component can flow without dissipation; and flow of the superfluid component is subject to 
quantum restrictions. One of the major such restrictions is that the vortices in superfiuids 
are quantized [6-11]. Emphasizing that the peculiarities of turbulence in superfiuids, as 
compared to that in classical fluids, are due to quantum effects, turbulence in superfiuids is 
called quantum turbulence [12-14]. 

There are several ways of producing turbulence in superfiuids, e.g., by realizing coun- 
terflows between the normal and superfluid components, by moving grids or by vibrating 
objects immersed in the superfluid [12-14]. The fully developed quantum turbulence is rep- 
resented by a random complex tangle of quantized vortices. Their appearance can be due to 
different reasons, to mutual friction between the normal and superfluid components [12-16], 
accompanied by the Kelvin-Helmholtz instability [17,18], to thermal and quantum tunneling 
[12], and to other dynamic instabilities [19,20]. When the external disturbance, producing 
turbulence, ceases to act on the fluid, turbulence decays through the Richardson cascades 
[12-14], the emission of sound radiation produced by vibrating vortices [21,22], by vortex re- 
connecion [23,24], and by Kelvin- wave cascades [25,26]. The evolution of quantized vortices 
can be visualized by injecting into the fluid small admixture particles [27]. 

Bose-Einstein condensed gases of trapped atoms [11, 28-33] provide additional possibil- 
ities for studying superfluids in a well controlled way. However, by rotating trapped Bose 
condensates, one can create quantized vortices, which tend to form crystalline structures, 
when the rotation velocity increases [34-36]. Simple rotation does not produce turbulence 
in trapped systems. In order to create a tangle of vortices, typical of quantum turbulence, 
it is necessary to use some other techniques. Actually, a vortex is a particular case of coher- 
ent topological modes, among which there exist many other soliton-like excitations [37,38]. 
Hence, the general question is how it would be possible to create coherent topological modes 
for trapped atomic systems? 

A method for generating various coherent topological modes in traps has been advanced, 
based on the modulation of atomic cloud density [39-41]. The modulation can be done 
in two ways, either by oscillating trapping potentials [39-43] or by alternating the atomic 
scattering length using Feshbach resonance techniques [33,44,45]]. Both these ways yield 
similar results, so that it is a matter of convenience which of them to employ. The generation 
of the coherent modes can be achieved by rather weak modulation amplitudes, but provided 
that the frequency of the alternating trapping potential or of the alternating scattering length 
is in resonance with the transition frequency corresponding to the desired coherent mode [39- 
45]. By increasing the modulation amplitude, it is feasible to realize the nonresonant mode 
generation [42] . In particular, the generation, without the direct resonance condition, can be 
done, e.g., by means of parametric conversion or harmonic generation [46,47]. Very strong 
forced oscillations can even destroy superfluidity at all [48]. The alternating modulation of 
trapping potentials, to some extent, can be modelled by combined rotations around several 
axes [49]. 

For theoretical description of quantum turbulence, it is sufficient to consider the nonlinear 
Schrodinger equation (Gross-Pitaevskii equation) [49,50], keeping in mind low temperatures 



and weak atomic interactions, when almost the whole system is in the Bose-condensed state. 
The generation of coherent topological modes and, hence, turbulence, can also be realized 
under more general conditions of finite temperatures and strong atomic interactions [51]. The 
presence of the normal component is not necessary for the appearance of quantized vortices. 
The latter arise because of the occurrence of a heteroclinic instability [38] in a nonequilibrium 
system, when the homogeneous atomic density becomes dynamically unstable with respect 
to its transformation to a nonuniform density comprising topological defects, like quantized 
vortices. This is analogous to the development of classical turbulence that is accompanied 
by a strongly nonuniform distribution of density [52-55] . In quantum turbulence, the nucle- 
ation of vortices, under the dynamic instability, is caused by developing unstable collective 
excitations, such, e.g., as quadrupole modes [19,20]. The typical feature of quantum turbu- 
lence is the random spatial distribution of quantized vortices [12-14], similar to the models 
of dynamical glass [56] . 

Experimentally, the modulation of the trapping potential has been used for the generation 
of quantized vortices [57] and vortex tangles [58]. The appearance of vortex tangles implies 
that the system enters the fully turbulent regime [12-14]. 

Turbulence is a strongly nonequilibrium phenomenon representing an essentially nonuni- 
form matter. This is why it is difficult for theoretical description. But, in the case of fully 
developed stationary turbulence, it can be possible to invoke statistical methods. It is the 
aim of the present paper to suggest a statistical model of stationary quantum turbulence. 

2 Vortex mixture 

Quantized vortices are the necessary ingredients of quantum turbulence. Each vortex can 
be represented as a filament passing through the fiuid. The velocity v = v(r) at a point r, 
due to a filament, is given by the Biot-Savart formula 

n f (s-r) 

X ds , (1) 



2m J |s — r^ 

where m is atomic mass, n is an integer characterizing the quantum of circulation, the 
vector s refers to a point on the filament, and the integration is taken along the filament. 
The filament has a definite direction associated with its vorticity 

a; = V X V . (2) 

Integrating velocity along a path around the vortex core gives the quantum of circulation, 

vrfl=^^ (n = 0,±l,±2,...) . (3) 

m 

For example, in the case of a rectilinear vortex, with the vorticity along the axis z, the 

velocity is 

n n 

V = e^ = — 2" e^ X r_L , 

mr± mr\_ 

where r±_ = \/x^ + -ip- = \ r±\. This gives the vorticity 

Un = 6{x)6{y)e^ . 

m 

3 



The case of n = means the absence of a vortex. 

In the turbulent superfluid there are many vortices with different vorticities, whose dis- 
tribution can be characterized by a measure m{u). Thus, in the case of rectihnear vortices, 
integration of a function /(cJ), with a measure m{u), gives 



f{u) dm{u) = ^ fi^n) dQ{e, 



where summation is over the circulation quanta n = 0, 1,2, . . . and integration is over the 
spherical angle around e^, that is, over all directions of vorticity. 

A fluid with a quantized vortex represents a dynamic state that is principally different 
from the state without vortices. Moreover, the states with the vortices, having different 
winding numbers or vorticities can be treated as different states. At any instant of time, 
the system, formed by multiple regions, each containing a vortex, can be considered as a 
composition of states corresponding to these different regions. So, at each instant of time, 
the system is represented as a configuration, composed of regions with different vortices. The 
density of vortices in the system and their directions are random in space and are varying 
in time. This means that the vortex configuration varies. The statistical description of such 
a situation assumes an averaged picture corresponding to the averaging over all admissible 
state configurations. The theory of describing such heterophase mixtures has been advanced 
in [59-63] and reviewed in [64-66], where all mathematical details of the related averaging 
procedure are elucidated. 

Following this approach [64-66], we come to an effective Hamiltonian characterizing the 
random mixture. In the present case, this Hamiltonian is the direct integral 



H= H{uj) dm{uj) H2 



(4) 



Here H{od) is the Hamiltonian of a system with a vortex characterized by the vorticity u 
and H2 is the Hamiltonian corresponding to the normal (nonsuperfluid) state. The latter 
has to be taken into account since, in the process of generating turbulence, superfluidity can 
be destroyed in a part of the system. This Hamiltonian (4) is defined on the mixture space 

M = (^n{u)(^n2, (5) 

which is a continuous tensor product of the weighted Hilbert spaces. More details on the 
meaning of continuous products are given in the Appendix. Each weighted space 'H(a;) is a 
copy of the system Hilbert space, with a weighted scalar product [64-66] , corresponding to 
the state with a vortex labelled by its vorticity u. And 7^2 is a space corresponding to the 
normal (nonsuperfluid) state. 

The superfluid states with the vortices are described by the grand Hamiltonians 

H{uj) = H{uj) - fioi'^)Noiuj) - /ii(cI;)iVi(a;) , (6) 

in which H{od) is the energy operator, Nq{uj) is the number of condensed atoms, A'^i(a;) is the 
number-of-particle operator for uncondensed atoms, and ^o{^) and /ii(cJ) are the Lagrange 
multipliers guaranteeing the normalization conditions for the number of Bose-condensed 



atoms No{u) and for the number of uncondensed atoms A^i(a;) in the superfluid state of 
vorticity u [51,67-69]. The nonsuperfluid state is characterized by the grand Hamiltonian 

H2 = H2- fi2N2 , (7) 

where H2 is the corresponding energy operator, N2 is the related number-of-particle operator, 
and ^2 is the chemical potential of the normal state. 

In what follows, we shall consider local atomic interactions 

$(r) = %5{r) , $0 = 47r ^ , (8) 

m 

where m is atomic mass and a^, scattering length. The local form (8) is not principal, but 
is taken just for the brevity of notations. Throughout the paper, the units are employed, 
where h = 1 and ks = 1. 

The Hamiltonian energy operator in (6) has the form 



H(ffj^wj(ffj / <i.t(ii,r) 



C + ^w 



■?/'(a;,r) dr + 



+ ^^^^i— ^ $0 / 'ip^{u,r)'ijj^{u,r)'ilj{u,r)'ijj{u,r) dr , (9) 



where the field operators ilj{u,r) act on the space 'H{uj) and the operator of momentum is 
p = —iV. The factor wi{u) is the geometric weight of the superfluid vortex state with 
vorticity u. The function f/(r) is an external potential, say, a trapping potential. The 
superfluid system is assumed to be Bose-condensed. And the Bose-Einstein condensation 
necessarily requires that the global gauge symmetry be broken [70]. 

The most convenient way of gauge symmetry breaking is by means of the Bogolubov shift 
[71], which can be represented as the sum 

tlj{uj, r) = r]{uj, r) exp{iS{r)} + ^i{uj, r) , (10) 

in which S{r) is a real function, such that 

V = v(r) = — VS(r) (11) 

m 

is the superfluid velocity. The first term in (10) is the condensate wave function and the 
second term is an operator of uncondensed atoms in the superfluid state with vorticity u. 
To exclude the double counting, the terms in sum (10) are to be orthogonal to each other, 

/ ■r]{u, r) exp{iS'(r)}-?/'|(a}, r) dr = . (12) 

And the statistical average for the operator of uncondensed atoms satisfles the property 

{M^,r)) = 0, (13) 

due to which the condensate wave function plays the role of the order parameter. 



The Hamiltonian energy operator in (7) is 



H2 = W2 ipl{r) 



l^ + ^w 



,2 



^ riS_ / J,Ur.\.l^l 



■?/'2(r) dr 



+ Y *o y ^Kr)^2(r)V^2(r)^2(r) dr , (14) 

where the field operators '?/'2(r), acting on 7/2, correspond to the normal (nonsuperfiuid) 
state and W2 is the geometric weight of the normal state. All field operators satisfy the Bose 
commutation relations. 

The condensate wave function defines the density of condensed atoms in the state with 
vorticity u, 

Po{uj,r) = \r]{uj,r)\'^ . (15) 

The density of uncondensed atoms in the superfiuid state, labelled by u, is 

pi{uj, r) = {ipliuj, r)^i(a;, r)) . (16) 

And the density of the normal phase is 

P2(r) = {4i^)M^)) . (17) 

The number of condensed atoms in the u vortex state is 

No{u) = w,{u) f po{u,r)dr. (18) 

The number-of-particle operator of uncondensed atoms in that state is 

Ni{u) = w^{u) j ij\{u,T)ij^{u,v) dr . (19) 

And the number-of-particle operator for atoms in the normal state is 

N2 = W2 j 4ir)Mr) dr . (20) 

So that the number of uncondensed atoms in the u vortex state reads as 

N,{u) = {N,{cu)) = w,iu) J p,iu,r) dr , (21) 

and the number of normal atoms being 

N2 = {N2)=W2fp2{r)dr. (22) 

The number of condensed atoms in all vortex states is the integral 

No= No{u) dm{u) , (23) 



while the number of uncondensed atoms in all vortex states is 

Ni= f Ni{uj) dm{uj) . (24) 

The total number of atoms in the system is 

N = No + Ni + N2. (25) 

The geometric weight of the superfiuid state with all types of vortices reads as 

wi = wi{uj) dm{uj) . (26) 



There can be in the system only two types of states, superfiuid and normal, which implies 
the normalization 

w,+W2 = l iO<w,<l), (27) 

where z/ = 1, 2. 

The equations of motion for the field variables are given by the related variations: for 
the condensate wave function of a superfiuid state with the vorticity u, 

for the field operator of uncondensed atoms in the superfiuid state with the vorticity u, 

and for the field operator of atoms in the normal (nonsuperfluid) state, 

^|:^2(r,t) = -^. (30) 

Let us introduce the notations for the anomalous average 

ai{uj, r) = {ipiiuj, r)^i(a;, r)) , (31) 

the anomalous triple correlator 

^{u, r) = (^l(a;, r)^i(a;, r)^i(a;, r)) , (32) 

and for the operator 

X{u, r) = 2tpl{u, r)tpi{uj, r)r]{u, r) + 

+ r]*{u, r)tpi{u, r)tpi{u, r) + tpl{u, r)tpi{u, r)tpi{u, r) , (33) 

where, for brevity, the time dependence is not explicitly shown. 

Then Eq. (28) yields the equation for the condensate wave function 



^ Q^V{^^^) = wi{u) 



'P + "^>'+f/(r)-Mc3) 



2m 



ri{uj,v) 



+ wl{uj)^o [poi'^, r)r]{uj, r) + 2pi{uj, r)r]{uj, r) + 0-1(0;, r)r]*{uj, r) + ^{uj, r)] . (34) 
Equation (29) results in the equation for the field operator of uncondensed atoms 






Wi[U 



^ + uir) - ,m 



V'i(^,r) + 



+ wl{u)^Q 2po(w,r)V'i(w,r) +77^(0;, r)?/'|(a;,r) +X(a;,r) 
And Eq. (30) gives the equation for the field operator of normal atoms 



' I ^^(^) 



W2 



2m 



+ f/(r) - ^2 



V^2(r) + W2'^oV^|(r)V^2(r)V^2(r) 



(35) 



(36) 



In these equations, again for brevity, the time dependence is not written explicitly. 

When considering a stationary state, it is possible to introduce an effective temperature 
Te/j characterizing the stationary input of energy into the system [72-74]. Then an averaged 
behaviour of the latter, in the sense of averaging over the observation time, can be described 
in terms of effective statistical ensembles and the corresponding thermodynamic quantities 
[75,76]. For stationary classical turbulence, a similar approach has been used by Kraichnan 
and Montgomery [77]. 

For the considered case of the stationary quantum turbulence, we can define the grand 
thermodynamic potential 

^ = -T^ff In Tr exp(-/3iJ) , (37) 

in which /3 = l/T^ff. The condition of stability for the heterovortex mixture requires that 
the thermodynamic potential (37) be minimal with respect to the corresponding geometric 
weights: 







> 



where, to take into account normalization (27), one has to set 

W2 = I — / uji{uj) dm{uj) . 
The first of conditions (38) gives 



(39) 



and the second yields 



6H 



5wi(uj] 



) = o 



''" }>Mi^)). 



6wl{u) 



6wi{uj) 



(40) 



(41) 



Since the right-hand side of inequality (41) is positive, the necessary stability condition reads 
as _ 

-^— > 0. 42 



Let us introduce the notations 



K(a)= iiia.i) 



l^ + ^W 



Tpiu^ r) ) dv 



$(c5) = $0 / ( 'ip\(^, 'c)i)\Q, r)-?/;(a;, r)-?/;(a;, r) ) dv 
for the superfluid state with vorticity u and 



(43) 



Ko 



4(r) 



£^+^« 



^2(r) ) dv , 



$2 = $0 / ( ^|(r)^|(r)^2(r)^2(r) ) rfr 
for the normal state. And let us define the integrals 



(44) 



Rq{Q) = / po(w, r) dv , -Ri(a;) = / pi(a;, r) lir , R2 = P2(r) (ir . (45) 

Then the numbers of condensed and uncondensed atoms in the u superfluid state, and atoms 
in the normal state are 

No{u) = wi{u)Ro{u) , Ni{uj) = Wi{uj)Ri{uj) , N2 = W2R2 , (46) 

respectively. Also, we define 

Ki{u) = K{u) - fio{u)Ro{u) - fii{u)Ri{u) , K2 = K2 - /i2i?2 • (47) 

Employing these notations, from Eq. (40), we find the geometric weight for the superfluid 
state with vorticity u in the form 



WiiUJ 



^2 + K2-Ki(uj) 



$(cJ) + $2 

/^From the stability condition (42) we obtain 

$(a;) + $2 > . 
The system chemical potential can be found from the relations 

F = n + ixN = n+[ fxo{uj)Noiuj) + fii{uj)Ni{uj) ] dm{uj) + fX2N2 , 

defining the free energy. This yields 



p 



fio{u)no{uj) + fii{uj)ni{uj) ] dm{uj) + /i2n2 , 



(4^ 



(49) 



(50) 



(51) 



where the corresponding atomic fractions are 






n2 



iV2 

N 



(52) 



The superfluid and normal phases are in mutual equilibrium, because of which their 
chemical potentials are connected with each other. The relation between the chemical po- 
tentials follows from the condition of equilibrium 5F = 0, where the variation is with respect 
to the numbers of atoms. This is equivalent to the equation 



5No{uj) + TTTT^ 5Ni{uj) 



5Nn(co 



5Ni{uj) 






5No 







Keeping the total number of atoms N fixed implies that 



[6No{u) + 6Ni{u)] dm{u) + 5N2 = Q 



By the meaning of the chemical potentials, 



5F 



6F 



Ho^uj) 



And for the atomic fractions, we have 



fXl[Uj) 



6N 



6N 



ni{uj) , 



SF 

Wo 



SN2 

6N 



/^2 



«2 



With the notation 



no = / no{u) dm{u) , rii = ni{uj) dm{u) 



normalization (25) becomes 

no + ni+n2 = I . 

Invoking the above relations, from the equilibrium condition (53), we obtain 

J[fio{uj)no{uj) + fii{uj)ni{uj)]dm{uj) 

f^2 : • 

no + ni 

Substituting this into (51) gives 

J[fio{uj)no{uj) + fii{uj)ni{uj)]dm{uj) 



/i 



no + ni 
that is. 

The total superfiuid density can be found [33] from the general expression 



Ps= P 



^ A^fPl 



3mV 



in which p = N/V is the total atomic density and the momentum dispersion 

A2(P) = (p2) - (p)2 



(53) 



(54) 



(55) 



(56) 

(57) 
(58) 

(59) 



(60) 
(61) 

(62) 
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characterizes the amount of heat dissipated in the system. The total momentum is defined 
as _ 

P = lim ^ , (63) 

where H^ is the Hamiltonian for the system boosted with velocity u (see details in [33]). 
This gives 

P = Pi + P2, Vi= rP{u)dm{u) , 



P(a;) = wi[u) I ip{od, r){—i'V)ip{uj, r) dr , 

P2 = W2 I Mr)i-i^)Mr) dr . (64) 

Because of the vanishing covariance 

COv(Pi, P2) = ^(PlP2 + P2P1) - (Pl)(P2) = , 

we have 

A2(P) = A2(Pi) + A2(P2). 

Therefore, the superfluid density (62) takes the form 

/3 



P2 = P - 



3mV 



A^(Pi) + A^(P2 



(65) 



Generally, a Bose system, subject to the action of an external alternating field, can be in 
the following five states, which depend on the modulation amplitude and time, that is, on 
the amount of energy pumped into the system. 

(i) Homogeneous superfluid: If there is no external perturbation, the system is in equilib- 
rium, and the effective temperature equals the usual temperature, Tg// = T. We assume that 
the latter is below the Bose-Einstein condensation temperature, so that the system is Bose 
condensed and is in a homogeneous superfluid state without vortices. The absence of vor- 
tices implies zero vorticity u = 0. The fact that the whole system is filled by a homogeneous 
superfiuid, without vortices, means that the geometric weight 

wi(0) = 1 . (66) 

This state remains under rather weak perturbations, when the effective temperature is yet 
close to T < Tc. 

(ii) Vortex superfluid: Disturbing the system by an alternating field injects energy into 
the system. The pumped mechanical energy increases the effective temperature. And, 
after sufficient amount of energy has been pumped into the system, there appear quantized 
vortices. First, there arise just a few of them, whose energy approximately equals the pumped 
energy. This can be expressed through the sum 

Wl{0) + Y,M^^) = ^, (67) 



11 



in which, in addition to the geometric weight wi{0), there appear a few terms wi{u), with 
cJ 7^ 0, describing a small (roughly speaking, less then ten) number of vortices. 

(iii) Turbulent superfluid: When the energy, pumped into the system, reaches a critical 
value, quantum turbulence develops with a great number of quantized vortices forming a 
tangle. This can be denoted as 

Wi{u) dm{uj) = 1 , (68) 

the number of vortices being much larger than ten. The integral includes the term ifi(O). 

(iv) Mixture of turbulent superfluid and normal fluid: A rather strong pumping, not only 
produces vortices, but starts destroying superfluidity, so that the admixture of the normal 
(nonsuperfiuid) liquid appears. This means that the equation 

wi{oj) dm{uj) + W2 = 1 (69) 

becomes valid, with a nonzero geometric weight of the normal fluid W2 > 0. 

(v) Normal fluid: When the pumping becomes very strong, it can destroy all Bose- 
Einstein condensate and, hence, superfluid, transferring the system into the normal state. If 
there exists solely the normal phase, then 

W2 = l. (70) 

The normal state can be turbulent, but this would be the classical turbulence. 
The overall picture is presented in the scheme of Fig. 1. 

3 Concluding remarks 

A statistical model is suggested, describing stationary states arising in a superfluid subject 
to the action of an external alternating field pumping energy into the system. The latter 
is represented as a continuous mixture of dynamic phases with different vorticities. The 
case of no vortices is included, corresponding to zero vorticity. The possible destruction of 
superfluidity by the pumping, leading to the developing normal (nonsuperfiuid) state is also 
taken into account. The vortices arise owing to dynamic instability corresponding to the 
appearance of unstable growing collective excitations. 

The first three regimes, shown in the scheme of Fig. 1, that is, homogeneous superfluid, 
vortex superfluid, and turbulent superfluid, have been realized, by modulating the trapping 
potential, in experiments [57,58] of the Bagnato group. A detailed analysis of these exper- 
iments and their relation to theory will be given in separate publications. Here, we give 
some estimates for clarifying the physical picture by presenting the typical parameters of a 
turbulent superfluid. 

In experiments [57,58], trapped atoms of ^^Rb are cooled down to form Bose- Einstein 
condensate. Hence, atomic mass is m = 1.445 x 10~^^ g and scattering length is a^ = 
0.577 X 10~^ cm. Almost all atoms have been condensed, composing the condensate of 
N = 2 X 10^ atoms, thermal fraction being rather small. A cylindrical harmonic trap is 
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characterized by the radial frequency uj^ = 2tx x 210 Hz = 1.319 x 10^ s ^ and longitudinal 
frequency cu^ = 27r x 23 Hz = 1.445 x 10^ s^^. The related oscillator lengths, 



^± = \l , h 



muj± V ^T^^. 



are l± = 0.744 x 10 ^ cm and Iz = 2.248 x 10 ^ cm. The average frequency and oscillator 
length, 



/ 2 \l/3 , /,2 7 \l/3 / ^ 

V TnojQ 



are uq = 0.631 x 10^ s^^ and /q = 1-076 x 10^^ cm. The effective volume of the condensate 

^e// = 7r/_L2/2 = 27r/o 

is Veff = 0.783 X 10^^^ cm^. Therefore, the average condensate density is p = N/Vefj = 
2.554 X 10^^ cm"'^. The average interatomic distance is a = p^^'^ = 0.732 x 10"^ cm. This 
shows that as <ti a and the gas parameters are small: pa^ = 0.491 x 10"'^ and p^^^ag = 0.079. 
That is, the atomic interactions are weak. 

The trap is subject to an external field modulation during the time text = 0.02 s —0.06 
s with an alternating potential of frequency Umod = 27r x 200 Hz = 1.257 x 10^ s~^. Thence 
the modulation period is tmod = '^'^/^mod = 5 x 10^^ s. The local-equilibrium time [29] 
tioc = m/hpag = 0.929 x 10^^ s is much shorter than the modulation period tmod, which 
implies that the system is always in local equilibrium. 

The typical size of a vortex core is the healing length 

1 1 



A/2mp$o y/SnpOs 

which is .^ = 0.519 x 10~^ cm, being of order of interatomic distance a. 
The sound velocity 



— $0 = — V47rpa, 
m m 

becomes c = 0.994 cm/s. This defines the coherence length l^oh = h/mc = 0.735 x 10^^ cm 
that is of order of the healing length. 

The vortex velocity (1), for n = 1, is of order 



2m/o \ ^ 

which gives v = 0.103 cm/s. 

Another typical velocity is that characterizing the velocity of atomic collisions Vcoi = 
h/mas = 12.65 cm/s during the collision time tcoi = ds/vcoi = 0.456 x 10^^ s. Thus, the 
relation between the characteristic velocities is 

t; < c < Vcoi ■ 

The characteristic vortex energy 
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is 5 = 1.912 X 10^ s-\ 

Summarizing the relations between the typical times yields 

And the relations between the characteristic lengths are 

Finally, the turbulent regime has to be formed by many vortices, N^or > 10. But there 
is the maximal number of vortices of order 

, X 2 ., X 2 

ly rt 



"max 



r^j 



I'coh J \ S 



before the vortices start essentially overlapping with each other. If this overlap becomes too 
strong, the superfluid state is destroyed. Therefore the number of vortices in a turbulent 
superfluid is in the interval 

10 < N.„„.. < 200 ^ 400 . 



't)or 



This is in agreement with experiment [58], where up to 200 vortices were observed. 

The turbulent superfluid, being a tangled mixture of many condensates with vortices, 
should exhibit a specific behaviour, when the trapping potential is switched off. Each par- 
ticular condensate, being released from the trap, expands anisotropically. But, since the 
mixture of many condensates, as a whole, is isotropic, all particular anisotropic directions 
are averaged out, so that the total turbulent cloud expands isotropically, keeping its aspect 
ratio during the whole free expansion. This effect was observed in experiment [58]. 

As is mentioned above, a more detailed analysis of these and related experiments will be 
given in separate publications. 

Acknowledgements. I am grateful to V.S. Bagnato for discussions of experiments [57,58]. 
Financial support of the Russian Foundation for Basic Research is appreciated. 
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Appendix. Continuous products 

The model of quantum turbulence in superfluids, suggested in this paper, is based on the 
notion of continuous products. For the self-consistency of the paper, in this Appendix, the 
main definitions are given, related to continuous products. Such products have been used 
for treating continuous heterophase mixtures [63,64] and models of continuous random walk 
[78,79]. 

Let a complex function /(x) : M — > C of a real variable a; G M be given. Suppose, 
we consider x in an interval [a, b] . Divide this interval on n parts of length Axi so that 
Y17=i ^^i = b — a. Define the limiting procedure 

n — 7> oo , Axi — > (i = 1, 2, . . . , n) . 

A definite continuous product of the function f{x) on an interval [a, b] is given by the 
limit 

b n 

l\f{x)^ \imT\f{a + Ax,), (Al) 



n— >oo 

1=1 



where the above limiting procedure is assumed. 

In order to transform this limit to a more tractable form, let us take the logarithm of the 
limiting equation, which gives 



log TT fix) = lim V" log f{a + Axi) . 



j=l 



The logarithm can be taken over any base. Notice that the right-hand side here is the 
definition of an integral, so that 



b rb 



log JJ/(a;) = / log/(a;) dm{x) , 

x=a •^"■ 



where m{x) is a measure on R. Since the logarithm can be taken over any base, we may 
chose the natural logarithm. Then exponentiation yields another form of the definition for 
a definite continuous product: 



I I f[x) = exp < / ln/(a;) dm{x) 



(A2) 



The logarithm of a complex function is given by the expression 

ln/(x) =ln|/(a;)| +i aigf\x) . 

Using the notion of indefinite integrals, it is possible to define an indefinite continuous 
product 



n-— »{/>■"" ")■ 



(A.3) 
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where c > is a positive constant. Then, we can introduce the action inverse to the 
continuous product as 



f{x) = exp 



1 m'(x) dx \ j 



where m'{x) = dm{x)/dx. 

It is straightforward to generahze the notion of a continuous product as follows. Let a 
complex function f{x) : D ^ C be given on a measurable domain D with a measure m{x). 
Then the continuous product of f\x) over the domain D is 






ln/(x) dm{x) \ . {AA) 



It is possible to introduce continuous products of spaces [63,64,80]. Let Hx^Zx G U^ be 
vectors of a unitary space Ux with a scalar product {ux-, Zx) and x G D. The manifold D is 
assumed to be measurable, with a measure m{x). The norm of Hx-, generated by the scalar 
product, is \\yx\\ = \/{yx,yx)- 

The continuous tensor product of unitary spaces, 

U = <^Ux (A5) 

xeO 

consists of the vectors that are continuous tensor products 

y = (^yx e U . (A6) 

x&i 



The scalar product in U is given by 



x&i ^-^^ J 

Respectively, the norm oi y eU, generated by the scalar product, is 



(A.7) 



When Ux is complete, it becomes a Hilbert space (complete unitary space). Then U is also 
a Hilbert space. 

Finally, one can introduce continuous tensor products of operators. Let Tix be a Hilbert 
space, with a; G D, where D is measurable with a measure m{x). And let an operator Ax be 
given on Tix- The operator norm can be defined as 

114 II II 

11^^11= sup — I 1— {y^^O). 

VxeHx Ill/a; 1 1 

The continuous tensor product of the Hilbert spaces Tix is 

"H = Q^ Tix , 

xen 
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similarly to (A. 5). An operator A on "H is an operator continuous tensor product 



whose action on y G "H is given by the tensor product 



Ayv = Q^ A,y, . (A.9) 

The operator norm of A G "H, defined in the standard way, 

Mil = sup M^ (y^O), 
yen \\y\\ 

results in the relation 

I |i| I =exp I Mn I |i^| I rfm(a;) I . (AlO) 

These definitions and properties of continuous products allow us to accomplish all neces- 
sary calculations for the considered turbulent superfiuid represented as a continuous vortex 
mixture, in which the role of the index x is played by the vorticity u. 
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i 

normal fluid 


■ ^'f^ w, = 1 


turbulent superfluid 
and normal fluid 


J w^(ub)dm(ub) + ^2 = 1 


turbulent superfluid 


J w^(ih)dm{ih) = 1 


vortex superfluid 


w,{0) + Xw{cv.)= 1 


homogeneous 
superfluid 


w,{0)=\ 



Figure 1: Scheme of the sequence of states for a superfluid subject to the action of an 
alternating external field. 
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